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Zaemit  report  considen  recent  experimental  findings  pertaining  to  the  vortex-excited  oscillations  of 
bluff  cylinders,  and  compares  various  approaches  to  measuring  the  resultant  fluid  forces.  As  one 
example,  measurements  of  the  structural  damping  and  croasflow  response  amplitudes  were  made  for 
several  systema  over  a lepietentative  range  of  wind  tunnel  flow  speeds.  The  fluid  dynamic  damping  or 
reaction  In-phaae  with  the  cylinder’s  velocity  also  was  meuured  for  different  Incident  flow  speeds 
over  a range  which  included  the  locking-on  of  the  vortex  shedding  to  the  resonant  cylinder  vibra- 
tions. A mathematical  model  It  propoted  in  connection  with  these  measurements,  and  the  components 
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0.  Abstract  (Continued) 

of  the  excitation  or  lift  force  on  a vibrating  structure  can  be  obtained  from  the  model  once  the 
structural  damping  and  fluid  reaction  forces  are  known. 

In  a second  approach,  the  fluid  forces  were  measured  on  a cylinder  which  was  forced  to 
oscillate  in  water  at  different  frequencies  and  crossflow  amplitudes  over  a range  of  Incident  flow 
speeds.  Both  the  steady  and  unsteady  forces  are  altered  significantly  in  the  locking-on  regime,  with 
large  amplifications  in  teth  the  lift  and  drag.  These  vortex-excited  forces  are  then  incorporated  into 
the  mathematical  model  for  self-excited,  resonant  vibrations  and  the  complementary  nature  of  the 
two  measurement  schemes  is  discussed.  The  results  obtained  from  the  present  study  are  suitable  in 
general  for  application  not  only  to  the  vortex-excited  vibrations  of  rigid  cylinders,  but  also  to  other 
flexible  bluff  structures  and  cables  in  air  and  in  water.  ^ 
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VORTEX-EXCITED  UNSTEADY  FORCES 
ON  RESONANTLY  VIBRATING,  BLUFF  STRUCTURES 

INTRODUCTION 

When  one  of  the  characteristic  frequencies  of  a bluff,  or 
unstreaml Ined , structure  is  near  the  frequency  at  which  vortices  are 
naturally  shed,  resonant  self-excited  vibrations  often  occur  if  the 
structure  is  lightly  damped.  There  is  also  a range  of  frequencies 
near  this  so-called  Strouhal  frequency  of  shedding  where  forced 
vibrations  of  a bluff  body  cause  the  vortex  frequency  to  be  captured 
by,  or  to  lock  onto,  the  body  frequency.  This  means  that  the  body  and 
wake  oscillations  have  the  same  frequency  and  that  the  Strouhal 
frequency  of  vortex  shedding  from  a stationary  body  is  suppressed. 

This  locklng-on,  or  wake  capture  as  it  is  often  called,  causes 
vortex-excited  oscillations  to  occur  over  a range  of  flow  speeds.  The 
steady  and  unsteady  forces  which  act  on  a structure  are  amplified  as 
a result  of  these  vibrations,  and  the  changes  in  these  forces  are 
closely  coupled  to  changes  in  the  flow  field  near  the  body. 

The  purpose  of  this  report  is  to  consider  recent  experimental 
findings  pertaining  to  the  vortex-excited  oscillations  of  bluff 
cylinders,  and  to  compare  various  approaches  to  measuring  the  fluid 
forces.  As  one  example,  measurements  of  the  structural  damping  and 
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crossflow  response  amplitudes  were  made  for  several  systems  over  a * 

i 
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representative  range  of  wind  tunnel  flow  speeds.  The  fluid  dynamic 
damping,  or  reaction  in-phase  with  the  cylinder's  velocity,  also  was 
measured  for  different  incident  flow  speeds  over  a range  which  included 
the  locking-on  of  the  vortex  shedding  to  the  resonant  cylinder  vibra- 
tions. A mathematical  model  is  proposed  in  connection  with  these 
measurements,  and  the  components  of  the  excitation  or  lift  force  on  a 
vibrating  structure  can  be  obtained  from  the  model  once  the  structural 
damping  and  fluid  reaction  forces  are  known. 

In  another  approach,  the  fluid  forces  were  measured  on  a cylinder 
which  was  forced  to  oscillate  in  water  at  different  frequencies  and 
crossflow  amplitudes  over  a range  of  incident  flow  speeds.  Both  the 
steady  and  unsteady  forces  are  altered  significantly  in  the  locking-on 
regime,  with  large  amplification  of  the  lift  and  drag.  These  vortex- 
induced  forces  are  then  incorporated  into  the  mathematical  model  for 
self-excited,  resonant  vibrations  and  the  complementary  nature  of  the 
two  measurement  methods  is  discussed.  The  results  obtained  from  the 
present  study  are  suitable  in  general  for  application  not  only  to  the 
vortex-excited  vibrations  of  rigid  cylinders,  but  also  to  other 
flexible  bluff  structures  and  cables  in  air  and  in  water. 

RELATED  INVESTIGATIONS 

A number  of  papers  and  reports  which  consider  various  aspects  of 
the  problem  of  vortex-excited  oscillations  and  their  forced-body 
analogues  have  recently  appeared,  as  have  the  proceedings  of  a recent 
symposium  on  flow-induced  structural  vibrations  [1].  The  topics 
discussed  at  this  symposium  ranged  from  basic  studies  of  vortex 
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formation  and  generation,  to  large-scale  field  studies  of  the  vibra- 
tions of  marine  structures,  and  to  architectural  aerodynamics.  For 
Instance,  Griffin  & Ramberg  have  examined  on  a relatively  small 
scale  the  effects  of  cylinder  oscillations  on  wake  similarity  f2], 
vortex  circulation  and  steady  drag  In  the  wakes  of  circular  cylinders 
at  normal  and  yawed  Incidence  [3,  4,  5].  King,  Prosser  & Johns  [6] 
have  reported  the  results  of  laboratory-scale  experiments  In  which 
they  Investigated  both  In-line  atid  crossflow  vibrations  of  flexible 
cylindrical  piles  In  steady  water  currents.  King  [7]  has  also  compared 
model  experiments  such  as  those  just  mentioned  with  full-scale 
experiments  on  offshore  structures  [8],  and  has  extended  his  Investi- 
gations to  Include  the  effects  of  Incident  flows  at  angles  other  than 
normal  to  the  axis  of  a flexible  cylinder  [9] . 

The  complexities  Involved  In  a complete  formulation  c 
equations  of  fluid  motion  past  stationary  and  vibrating  bluff  bodies 
still  preclude  a solution  from  first  principles,  and  so  alternative 
and  less  complicated  methods  have  been  sought  for  engineering 
applications.  A formulation  which  couples  an  equation  for  the 
structural  motion  with  a nonlinear  fluid  oscillator  for  the  fluctuating 
lift  was  first  proposed  by  Bishop  & Hassan  f 10]  and  formulated  with  a 
van  der  Pol  oscillator  for  the  lift  by  Hartlen  & Currie  [11].  Bishop 
& Hassan  were  the  first  to  observe  that  the  process  of  locklng-on 
involves  rapid  phase  shifting  of  the  unsteady  fluid  forces  relative  to 
the  vibration  displacement,  and  that  this  phase  transformation  is 
accompanied  by  large  amplifications  of  both  the  lift  and  drag  forces. 
Soon  after,  a modified  form  of  this  original  concept  which  more 
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accurately  predicted  the  resonant  response  of  an  elastically-mounted 
cylinder  was  proposed  and  compared  with  experiments  by  Skop  6e 
Griffin  fl2]  and  Griffin,  Skop  & Koopmann  [13].  Iwan  & Blevins  [14] 
also  have  proposed  an  oscillator  formulation  which  is  based  upon  the 
Integral  momentum  equation  for  a system  composed  of  a bluff  cylinder 
In  a uniform,  steady  flow.  These  so-called  non-linear  or  "wake- 
oscillator"  models  have  been  generalized  by  Skop  6e  Griffin  [15], 

Blevins  & Burton  [16]  and  Iwan  [17]  to  include  the  vortex-excited 
oscillations  of  elastic  structural  elements.  The  formulations  of 
Skop  & Griffin  and  Iwan  were  arrived  at  independently  but  are  identical 
and  yield  the  same  results  for  a variety  of  cylindrical  structures 
such  as  flexible  cantilevers,  taut  cables  and  pivoted, rigid  members 
in  air  and  water. 

A wake-oscillator  type  of  formulation  recently  was  applied  by 
Blake  [18]  to  the  problem  of  predicting  the  vortex-excited  response 
of  blunt-based  hydrofoil  struts.  Good  agreement  was  obtained  for  the 
predicted  and  measured  modal  displacement  amplitudes  of  a cantilever 
strut  under  varying  structural  damping  and  flow  conditions.  This 
successful  extension  of  the  nonlinear  oscillator  concept  to  the  case 
of  the  vortex-excited  bending  vibrations  of  hydrofoils  gives  further 
evidence  of  the  potential  importance  of  such  methods  in  engineering 
practice. 

Most  recently  the  unsteady  lift  forces  on  vibrating  cylinders 
were  measured  by  Griffin  & Koopmann  [19]  in  a wind  tunnel,  and  by 
Mercier  [20]  and  Sarpkaya  [21]  in  water  tunnels.  Further  evidence 
was  given  in  these  experiments  of  the  highly  nonlinear  dependence 


between  the  various  fluid  dynamic  force  components  and  the  amplitude 
and  frequency  of  the  vortex-excited  and  forced,  resonant  vibrations. 


UNSTEADY  FORCES  ON  RESONANTLY  VIBRATING,  BLUFF  STRUCTURES 

The  equation  of  motion  for  a resonantly  vibrating,  cylindrical 
structure,  excited  by  vortex  shedding,  can  be  written  in  the  form 


^ 2 
where  ( •)  = — . Here  y » y/D,  t = u)  t , u = and  is  the 

M 

structural  damping  ratio.  The  force  coefficients  are 


Lift:  C^  » slnioT  + rp),  o = ^ , (2a) 

^ PV  D n 


Reaction:  C » = — 

^ \ PV^D 


C ^ sin  (OTT  + cpj) , o = — , 

n 


as  defined  by  Griffin  and  Koopmann  fl9],  where  cp  and  cp^  are  the  phase 

angles  between  the  lift  and  the  displacement  and  the  reaction  and  the 

acceleration,  respectively.  The  two  forces  represented  by  C and  C 

U K 

are  orthogonal,  as  will  be  shown  later.  An  analogous  formulation  of 


the  equation  of  motion  recently  has  been  introduced  by  Blake,  Maga  & 
Finkelstein  (22],  who  investigated  the  vortex-excited  oscillations  of 
propel lore  and  hydrofoils. 


If  the  displacement  y is  written  as 


y - Y sin  OfT,  » = ~ 


then  the  equation  of  motion  separates  into 


sinCKT;  - a Y + Y-  p, 


(^)  [' 


C cos  cp  - C cos 
L R 


"ij- " 


cos  QfT;  2 C Y 

S 


sin  cp  - C si 
L R 


n cp^  j =0 


when  the  coefficients  of  sin  or  and  cos  ar  are  grouped  appropriately. 
The  various  force  components  are  Identified  as  follows 
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- O'  Y + 


STIFFNESS 
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ADDED  MASS 


/111  r 

k)  t 


COS  cp  - 


C_  cos  cp, 


STRUCTURAL 

DAMPING 


+ 2 C Y 
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FLUID 

EXCITATION  DAMPING 


C sin  cp 

L/ 


Cr  sin 


Using  a different  notation  Sarpkaya  [21]  has  expressed  the 
measured  fluid  force  on  a vibrating  cylinder  as  the  sum  of  orthogonal 
components  in  the  form 


SotAL  ■ * ^mh  - ^dh  ^ ^5) 

- pV  D 

where  C . is  an  "inertia"  force  and  C.,  is  a "drag"  force.  These 
mh  dn 


force  components  are  related  to  Sarpkaya's  "generalized  Morison 
coefficients"  and  by 


dh 


3n  ^di.  y b j r 


and 


('6a) 


(fib) 


2rv 

where  V is  the  "reduced  velocity"  V = — . 

r r (J)  D 


In  this  case  the  equation  of  motion  for  the  cylinder  is  written  as 


y + 2 Cg  y + y = u(  — 


■{if  !■ 


, , sin  CKT  - C„  cos  CVT  . (7) 
mh  dh 


1 


iThe  drag  force  C,,  is  negative  when  energy  is  transferred  to  the 
dh 

cylinder,  as  is  the  case  for  resonant,  vortex-excited  oscillations.) 
If  a steady-state  response  is  assumed,  then 

y * Y sin  fo'  t - f ) , 

where  C is  an  undefined  phase  angle.  This  equation  of  motion  then 
separates  into 

2 

sin  O'  t:  O'  Y cos  f + 2 C Y sin  f + Y cos  f 

s 


(y 


"mh  = " 


fSa) 


and 


COB  O'  t:  O'  y sin  ( + 2 C Y cos  c - Y sin  ( 
^ s 


•(<)' 


"dh-° 


(8b) 


when  the  coefficients  of  sin  a r and  cos  a t are  again  grouped  appro- 


priately. If  the  first  equation  is  multiplied  by  sin  € and  the  second 
by  cos  € , then 


2 a C Y - u 
s 


fi)’  t 


C . sin  C - C,-  cos  C 
mh  dn 


]■“ 


(9a) 


when  the  two  results  are  added.  Next,  equation  (8a)  is  multiplied  by 
cos  € and  equation  (8b)  by  sin  C,  to  obtain 


- O'  Y + Y - p 


, V 2 r 

/ u)^  \ r 

1 0)  y 1 * 


C , cos  € + C,,  sin  € 
mh  dh 


0 (9b) 


when  the  two  equations  are  subtracted.  The  terms  in  equations  (9a) 
and  (9b)  can  be  arranged  and  labeled  in  the  same  manner  as  before: 


STRUCTURAL 

INERTIA  STIFFNESS 

- Y + Y » 


(</  t 
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C..  sin  e + C , cos  € 
dh  mh 


STRUCTURAL 

DAMPING 
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ii)  [- 


FLUID 
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C..  cos  C + C , sin  ( 
dh  mh 


These  equations  have  the  same  form  as  those  formulated  by  Griffin  & 
Koopmann  [19],  and  the  two  approaches  are  in  fact  identical  when 


sin 

cp  . 

-Cdh  cos  e. 

excitation'*' 

(10a) 

c 

L 

cos 

cp  > 

^dh 

FLUID  INERTIA 

(10b) 

^R 

sin 

'Pi 

- -‘^mh 

FLUID  REACTION  (DAMPING) 

(10c) 

^R 

cos 

'Pi 

ADDED  MASS 

(lOd) 

'^See  Appendix  I for  a further  discussion  of  these  equations. 
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These  relations  can  be  used  to  compare  results  from  the  two  approaches. 


From  the  above  system  of  equations  It  follows  that  tan  cp^  = tan  f 


tan  cp 


1 

tan  cp, 


(Recall  that  € and  cp  are  measured  from  y,  whereas  cp^  Is  measured 
from  y.)  It  also  follows  from  these  equations  that 


and  that 


2 2 


L'  ’ dh' 


• - I^rI  - 

The  phase  relations  between  the  various  forces  components  are  shown  In 
Figure  5(b)  for  a typical  set  of  conditions. 

A relation  between  the  displacement  Y and  the  total  fluid  force 
can  be  obtained  by  squaring  and  adding  equations  (4a)  and  (4b)  to  yield 


(l-c.2)  y2^^2  /^\  + 


- 2 (cos  cp  cos  cp^  + sin 


ft) 

+ sin  cp  sin  cp^) 


This  equation  simplifies  to 


1^(1  - cy)2+  (2C30')2jY2  ’ [ 


C 2 2 

L R 


since  tan  cp  ■ - , The  phase  angles  cp  and  cp,  can  be  retrieved 

tan  cp^  ® 1 

from  the  equations 


2 C Y 
s 


^T  r r 'P  ■ IT  “In  , 

\ n / L T T J 


r 


(1  - a ) Y - p. 


. 2 
\ 

C? 

— cos  cp 
Ct 

- ^ cos  CPj^ 

where 


Also  choose  - cp  = n/2,  so  that  these  equations  reduce  to 

()2 

r r-  sin  p cos  cp,  - cos  P sin  cp J 

1 


2 a C V = p. 

s 


and  2 

(1  - ^2)  Y = [+  sin  P sin  cp^  - cos  P cos  cp^l 


where  P = arc  tan  ^ . Simplifying  further,  one  obtains  the  equations 


/ U)  \ ^ 

- ^ “ 2 a Cg  Y, 

- p,  C,j.cos(cp^  + p)  = (1  - » ) Y, 


and 


tan  (cp.  + P) 


2 ^gs 

(1  - Ot^) 


W. 


(13a) 


(13b) 


(13c) 


This  latter  equation  can  be  rearranged  as 


W - tan  P 

tan  cPi  - j ^ p 


(13d) 


If  W - tan  Y*  then 

tan  cpj^  ■ tan  (y  ■ P) 


(14a) 
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and 


I 


cPj  = V - P. 


The  tangents  of  the  various  phase  angles  are 


(14b) 


'^R  “ ‘^L 


(15a) 


tan  cp  - 


1 

tan  cp 


- W 


(15b) 


or,  alternately,  the  angles  cp  and  cp^  are 


cp  - arc  tan 


cp^  - arc  tan 


^ V 

\ “ <^R  / 

(Cl  + « C,  \ 


(16a) 


(16b) 


according  to  the  formulation  of  Griffin  & Koopmann  [1]. 

Similar  expressions  can  be  found  from  the  formulation  employed 
by  Sarpkaya  [21]  and  his  students  {see  Raposo  [23],  Meyers  [241  and 
Fortlk  [25l}.  The  equations  of  motion  reduce  to 


ik)  [s 

> - V . R (^)  c,  [; 


C u C.. 

mh  , - dh 
- — sin  € - - — cos 

Ct 


cos  € + 


^dh  sin  € 

"t  J 


where 


'c  2 + c 2 

^mh  dh 


'1 


“iiih 


If  sin  and  cos  - 7;^.  then  these  equations  become 


(1  - ry  ) Y - p, 
which  yield  the  result 


ft)' 

(J) 


sin  (€  - P ) 


cos  (C  - P ') 


2 O'  C Y 2 a C 
a s s 

tan  (f  - P ) = ~ 


(I  - cy2)  Y (1  - 0-2) 

This  equation  can  be  rearranged  as 


W. 


(17a) 


(I7b) 


(17c) 


tan 


, W + tan  P ' 


1 - W tan  p 


(18) 


If  W = tan  V as  before,  then 


where 


and 


f ■ V + P ^ . 

c 


P » arc  tan 


Y “ arc  tan 


(C) 


(19a) 


(19b) 


In  order  for  self-excited  oscillations  to  take  place,  energy  must 
transferred  to  the  cylinder  and  Its  mounting  system  from  the  fluid. 

The  energy  transfer  over  a cycle  of  the  oscillation  Is  given  by  the 
Integral 


■/ 


y*  Fl  y dt 


(20) 


The  lift  force  Is,  In  Sarpkaya's  notation. 


\ (1/2  P V D L ) cos  or  T , 


+ See  Appendix  1 for  a further  discussion  of  the  work  done  and  energy 
dissipated  by  the  various  forces  over  a cycle  of  the  cylinder  motion. 
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and  the  cylinder  vibration  velocity  is 


y = O'  Y D cos  (O'  T - e)  , 
so  that  the  equation  for  the  energy  transfer  is 


- 1 2 2 
E - i p D L 


U)  T 

f " ^‘^dh^ 

Jo 


(Y  O')  cos  (or  T - €)  dT 


2 2 

1 p V D L 


Y I cos  u cos  (u  - €)  du, 


n (-C..  ) Y cos  e , 
an 


(21a) 


A necessary  condition  for  self-excited  oscillations  to  take  place  is 
E . P <-Cj^)  Y CO.  t > 0 

or,  equivalently,  cos  € > 0 since  is  positive.  This  leads  to 

the  condition 


C u - W C ,,  > 0 
mh  dh 


(21b) 


for  self-excited  oscillations  since 


cos  € - 


<=0^  - » ‘^dh 
- ” '=ch>'  + «=dh  * “ 


In  the  notation  adopted  by  Griffin  & Koopmann  (191,  the  energy 
transfer  is 


E - 1/2  p L 


/.U)  1 

f" 

•'A 


sin  (o  T + cp)  (Yo)  cos  o/r  dr 


E-C^Y 


sin  (u  + «p)  cos  u du  ■ n C,  Y sin  V.  (22) 


i ^ 

f i-' 

1 

' , , .4,13  -i 

! 1 t J 

-W‘ 

I " / ’ 


Since  Is  positive,  a necessary  condition  for  self-excited  oscilla- 
tions is  sin  > 0.  This  leads  to  the  condition 

C + W C,  >0 
R L> 

since 

sin  cp  = — ■ ' - - " - ' ' ' ' • 

* <‘l  - "V 

This  is  analogous  to  the  condition  derived  between  and 

UNSTEADY  FLUID  FORCE  MEASUREMENTS 

Some  interesting  and  useful  comparisons  can  be  made  between  the 
fluid  forces  measured  when  a cylinder  is  forced  to  vibrate  and  the 
fluid  forces  measured  when  the  cylinder  is  resonantly  excited  by 
vortex  shedding.  The  forced  cylinder  measurements  discussed  here 
were  made  by  Sarpkaya  [21]  and  Mcrcler  [20],  and  the  self-excited 
cylinder  measurements  were  made  by  Griffin  & Koopmann  [19]  and  others 
[26,  27,  28] . 

Typical  examples  of  the  measurements  reported  by  Sarpkaya  are 
shown  in  Figures  1 and  2.  The  measured  values  for  the  inertia 
coefficient  are  shown  in  Figure  1 as  a function  of  the  reduced 
velocity  V^,  for  a displacement  amplitude  of  Y - 0,50.  Of  particular 
note  is  the  rapid  change  in  in  the  vicinity  of  ■ 5,  This 
effect  corresponds  to  the  onset  of  locking-on  between  the  vortex  and 
vibration  frequencies  when  the  characteristic  frequency  of  the  flow 
approaches  the  vibration  frequency.  The  drag  or  resistance  force 
is  plotted  in  Figure  2 as  a function  of  at  the  same  displace- 
ment amplitude,  Y ■ 0.50.  This  component  of  the  total  fluid  dynamic 
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force  becomes  negative  near  « 5.  The  terminology  employed  by 
Sarpkaya  Is  somewhat  unwieldy  and  Is  an  extension  of  his  previous 
work  [29]  with  the  so-called  Morison  equation  [30]  for  purely  periodic 
flow  past  stationary  cylinders.  In  effect,  is  the  negative  of  the 
lift  coefficient  as  it  is  usually  characterized,  so  that  the  negative 
value  of  near  = 5 suggests  a net  transfer  of  energy  to  the 
cylinder  in  that  region.  These  forced-cylinder  results  can  be 
compared  to  the  vortex-excited,  transverse  oscillations  of  structures 
which  typically  occur  in  the  range  = 4.5  to  8 when  the  damping  of 
a structure  is  sufficiently  small  [19,  26,  31]. 

The  excitation  component  of  the  lift  force  C,  sin<>  = -C.,  cos  f 

L an 

measured  by  Griffin  & Koopmann  [19]  is  plotted  in  Figure  3 for  two 
cases.  For  the  first  case,  termed  System  I,  the  cylinder  reached 
a maximum  displacement  amplitude  of  Y • 0.47,  while  for  the  second 
case  the  maximum  displacement  was  Y <•  0.27.  The  conditions  for  these 
experiments  are  listed  in  Table  1.  The  lift  force  coefficients  in 
Figure  3 were  deduced  from  structural  and  fluid  dynamic  damping  or 
reaction  forces  on  freely-vibrating  cylinders  in  a wind  tunnel.  A 
typical  example  of  the  fluid  dynamic  damping  ratio  Cr-  ^ function 
of  the  reduced  velocity  is  shown  In  Figure  4.  There  is  a marked 
decrease  in  the  damping  in  the  region  of  the  locking-on,  and  much 
the  same  effect  has  been  observed  by  Blake,  Maga  & Flnkelstein  [22] 
during  measurements  of  the  vortex-excited  bending  vibrations  of 
blunt-based  hydrofoil  models. 

A comparison  now  will  be  made  between  the  force  components  of 
Sarpkaya' s forced  vibration  experiments  and  the  vortex-excited  forces. 
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Consider  as  an  example  the  measurements  » -0.5  and  » -0.2 

at  » 5.5  and  Y « 0.25,  reported  by  Sarpkaya  [21],  and  assume  W = 0.05 

as  a typical  value  (see  reference  19].  For  this  case  f » 204.7°  and 

C u - W C - -0.5  - 0.05  (-0.2)  - - .49, 
mh  dn 

so  that  self-excited  oscillations  do  not  take  place  (see  equation  21b). 

A phase  diagram  of  the  various  force  components  is  shown  as  Case  1 In 

Figure  5(a).  As  an  additional  example  consider  a set  of  conditions 

such  that  C^j^  » 1.6  and  C^j^  » -0.9  at  = 5 and  Y = 0.5.  Assuming 

once  again  for  convenience  that  W * 0.05,  then  f = -26.5°  and 

C , -W  C^^  - 1.6  - 0.05  (-0.9)  = 1.64, 
mn  an 

so  that  self-excited  oscillations  are  possible  for  these  conditions. 

This  Is  shown  as  Case  2 In  Figure  5(b) . 

The  excitation  component  of  the  lift  force  Is  defined  as 

c'  ■ C sln<^  ■ -C,,  cos  e (10a) 

L L an 

and  Is  Important  from  the  standpoint  of  energy  transfer  to  the 

vibrating  structure.  The  measurements  made  by  Sarpkaya  can  be  compared 

to  the  measurements  of  other  Investigators  by  means  of  the  above 

equation.  The  maximum  value  of  the  force  coefficient  in  Figure  2 

occurs  at  - 5,  and  the  same  result  Is  obtained  at  the  displacement 

amplitudes  Y > 0,13,  0,25  and  0.75  [21],  The  several  values  of 

-C,.  and  C , are  listed  In  Table  2 together  with  the  related  values  of 
dh  mh 

the  force  components  cos  C,  For  all  of  these  cases  of  forced 


vibration  the  condition  for  self-excitation 


excited  oscillation.  The  results  for  -C  cos  € are  plotted  against 

dn 

the  peak  displacement  In  Figure  6 together  with  a host  of  similar 
findings  from  the  results  obtained  by  other  Investigators.  The 
conditions  under  which  the  experiments  were  performed  are  described 
In  Table  3. 

Mercler  [20]  measured  the  force  coefficients  C . and  C..  on  a 

mh  dh 

cylinder  which  was  forced  to  vibrate  In  water.  These  measurements  were 

not  nearly  as  extensive  and  precise  as  were  those  of  Sarpkaya,  but  it 

Is  possible  to  estimate  the  various  force  components  given  by  equations 

(lOa  to  c)  from  Herder's  published  data.  The  results  are  listed  In 

Table  2(b)  and  are  plotted  In  Figures  6 and  7.  The  minimum  value  of 

C.,  was  measured  once  again  at  V = 5.5-6,  and  for  all  of  the  cases 
an  r 

listed  In  Table  2(b)  the  condition 


C . - W > 0 
mh  dh 

was  satisfied.  The  estimated  force  coefficients  -C..  cos  f from 

dh 

Herder's  measurements  are  In  reasonably  good  agreement  with  the  other 
results  plotted  In  Figure  6.  All  of  the  lift  force  measurements  in 
the  figure  exhibit  the  nonlinear  dependence  upon  amplitude  that  Is 
characteristic  of  vortex-excited  oscillations.  This  nonlinear 
"wake-oscillator"  was  first  observed  by  Bishop  & Hassan  f 10] , and  has 
since  been  studied  in  considerable  detail  by  others  [11,  12,  14]. 

The  measured  vortex-excited  displacements  In  Figure  6 are  scaled 


by  the  modal  response  factor 


which  takes  account  of  the  normal  modes  t^(x)  of  a flexible  structure 
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of  Inonersed  length  L.  As  shown  by  Skop  & Griffin  [151  ®nd  Iwan  [17], 


I 


1 


(x)  dx 


(x)  dx 


(23) 


For  the  special  case  of  a circular  cylinder  l'l'j^(x)|  » » 1,  and  the 

corresponding  values  for  cantilevers  and  pivoted  rods  are  listed  in 
reference  31. 

Several  important  characteristics  of  the  lift  force  behavior 
during  locklng-on  are  clear  from  the  results  plotted  in  Figure  6. 

First,  there  is  a maximum  of  the  excitation  component  of  the  lift  force 
at  a peak-to-peak  displacement  between  0.6  and  1 diameters  for  all  of 
the  cases  shown  in  the  figure.  Second,  the  maximum  of  the  lift 
coefficient  is  approximately  0.5  to  0.6  for  most  cases,  a notable 

exception  being  Sarpkaya's  result,  c'  ■»  0.75.  The  excitation  force 

Lj 

coefficient  decreases  to  zero  at  displacement  amplitudes  near  Y * 1, 
and  this  results  in  a limiting  amplitude  of  vortex-excited  oscillation 
of  one  to  one  and  one-half  diameters.  Sarpkaya's  measurements  give 
convincing  evidence  of  this  limit,  since  in  his  water  tunnel  experi- 
ments with  forced  cylinders  the  sign  of  the  force  coefficient 
changed  from  negative  to  positive  as  the  vibration  displacement  was 
Increased  from  0.75  diameters  to  1.03  diameters.  Mercier  [20]  also 

observed  a change  in  the  sign  of  the  force  coefficient  C..  as  the 

on 

vibration  displacement  y was  increased  to  1.3  diameters  from  the  values 
listed  in  Table  2(b) . 

The  formulations  for  the  fluid  dynamic  forces  due  to  vortex 
shedding  introduced  by  Blake,  Maga  & Finkelstein  [22]  and  Griffin  & 


i 
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Koopmann  fl9]  are  baaed  on  the  hypothesis  that  both  flow-induced  lift 
(excitation)  and  reaction  (damping)  forces  act  on  a resonantly 
vibrating,  bluff  structure.  The  fluid  damping  forces  were  measured 
In  air  by  Griffin  & Koopmann,  see  Figure  U,  and  In  water  by  Blake  and 
his  colleagues,  and  the  lift  coefficient  was  obtained  by  equating  the 
lift  component  C sin  cp  to  the  sum  of  the  structural  damping  and  fluid 
reaction  forces  as  given  by  equation  (lOa). 

It  Is  possible  to  deduce  the  damping  or  reaction  effect  of  the 
fluid  from  Sarpkaya's  total  force  coefficients  by  means  of  equation 
(10c) , 

^1  ■ '^mh  "• 

as  a further  step  in  comparing  the  two  approaches  described  in  the 

previous  section.  The  results  are  listed  in  Table  2(a)  and  are 

plotted  In  Figure  7,  and  they  yield  further  evidence  that  the  fluid 

force  decomposition  described  by  equations  (10a)  through  (lOd)  is  an 

appropriate  one.  A comparison  between  the  force  components  deduced 

from  Sarpkaya's  forced-vibration  measurements  of  C , and  C^,  and 

mh  dh 

plotted  In  Figures  6 and  7 shows  that  the  excitation  and  reaction 

forces  are  very  nearly  equal.  Though  Herder's  force  measurements 

were  not  as  extensive  as  those  of  Sarpkaya,  it  Is  possible  to  estimate 

the  reaction  effect  of  the  fluid  from  his  faired  curves  of  C . and  C ..  . 

mh  dh 

The  results  are  listed  in  Table  2(b)  and  plotted  In  Figure  7.  The 
same  reaction  or  damping  forces  component  as  reported  by  Griffin  & 

Koopmann  Is  plotted  in  Figure  7,  and  all  of  the  measurements  show  the 
same  general  pattern  of  behavior  even  though  these  latter  experiments 
were  performed  In  a wind  tunnel. 

I 


It  is  important  to  remark  upon  the  complementary  nature  of  the 
two  approaches  just  described  for  measuring  the  vortex-excited  forces 
on  bluff  bodies.  For  instance,  the  Morison-type  force  coefficients 


C . and  C , as  obtained  from  a Fourier  decomposition  of  the  total 
mh  dh 

force  measured  by  driving  the  cylinder,  do  not  give  the  necessary 
detailed  picture  of  the  various  force  components  which  act  due  to 
resonant,  vortex-excited  oscillations.  On  the  other  hand,  it  is  an 
exceedingly  difficult  experimental  undertaking  to  measure  the  resultant 
fluid  dynamic  force  on  a resonantly  vibrating,  self-excited  bluff 
cylinder.  Only  the  fluid  damping,  structural  damping  and  added  mass 
components  of  the  total  force  system  are  conveniently  measured  in  this 
latter  case.  Thus,  in  applying  the  results  from  the  two  approaches 
just  discussed,  it  is  necessary  to  recognize  the  limitations  Inherent 
in  each  and  to  take  advantage  of  their  complementary  nature  in  solving 
problems  associated  with  the  vortex-excited  vibrations  of  bluff 
structures . 
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L. 


It  is  useful  to  consider  the  work  done  and  the  energy  dissipated 
over  a cycle  of  the  cylinder's  oscillation  by  the  various  forces  in 
equations  (1)  and  (7).  The  work  and  energy  dissipation  are  given 
by  the  well-known  equation 


AE 


- f"  - - c 

W = / F dy  = I 

o •'o 


dt 


dt. 


(Al) 


Employing  the  nomenclature  Introduced  earlier,  and  evaluating  this 
Integral  for  the  structural  and  fluid  dynamic  forces  in  equation  (1), 
one  obtains 


2 C, 


UD  T 

/■ 


(y)^  dT  = M. 


(<) 


/•UU  T 

/•: 


sin  (o  T + cp)  y dT 


/•u)  T 


Cj^  sin  (O'  T + y di 


The  structural  inertia  and  stiffness  terms  are  conservative  and  do 
no  net  work  over  a cycle  of  the  oscillation.  The  structural 
velocity  y is 


y - O'  Y cos  (o  t)  , 
so  that  equation  (Al)  reduces  to 


or 


(y)  dx 


•Ol-^ 


Y gin  cp  - TT  Co  Y sin  cp 
K 1 


When  equation  (7)  Is  evaluated  In  the  same  manner,  one  obtains 


2 Cs  1 


I T ru)  1 

jj  ^mh  ^ y dx 


cos  O'  X y dx 


From  the  assumed  solution  to  equation  (7), 


y » o Y cos  (o  X - €) 


so  that  equation  (A3)  reduces  to 


2C./  " 

o V n ' 


C Y I 

mb  / sin  u cos  (u  - C)du 


- Y y cos  u cos  (u  - €)  du 

o 


c.  r 

o ' n ' 


n C . Y sin  C 
mh 


- n C Y coa  € 
an 


A straight-forward  comparison  of  equation  (AA)  with  equation  (A2) 
yields  the  results 


C sin  cp  = - cos  € EXCITATION  (10a) 

L dh 

C„  sin  cp  - - C , sin  C FLUID  REACTION  (DAMPING)  (lOc) 

R 1 mh 

as  previously  had  been  obtained.  The  "generalized  Morlson  coefficient" 

C , from  Sarpkaya's  formulation  [21],  has  a component  C , sin  C which 
mh  mh 

acts  to  oppose  the  velocity  y of  the  structure  (since  c Is  typically 

negative  for  vortex-excited,  resonant  oscillations)  and  acts  effectively 

as  a fluid  reaction  force  as  had  been  hypothesized  by  Griffin  & 

Koormiann  [19].  Likewise,  the  Morlson  coefficient  C , which  Is 

dh 

negative  during  vortex-excited  oscillations,  as  shown  In  Figure  2, 

is  In  effect  a lift  force  which  transfers  energy  to  the  vibrating 

structure  through  the  component  cos  €.  This  latter  force  Is  in 

phase  with  the  velocity  term  y.  The  results  given  by  equations  (A2) 

and  (AA)  show  that  the  excitation  force  on  the  cylinder  Is  balanced 

by  a fluid  reaction  force  and  a str’jctural  damping  force.  The 

reaction  or  damping  effect  of  the  fluid  can  be  measured  simply  In 

terms  of  an  effective  fluid  damping  factor  Cr.  shown  In  Figure  A, 

r 

which  is  taken  from  the  wind  tunnel  experiments  with  resonantly 
vibrating  circular  cylinders  reported  by  Griffin  & Koopmann  [19]. 

Blake,  Maga  & Flnkelsteln  [22]  have  made  analogous  measurements  In  a 
water  tunnel  with  model  hydrofoils  which  were  resonantly  excited 
by  vortex  shedding. 
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Figure  2 

plotted  against  V for  Y = 0.50;  from  Sarpkaya  [21]. 


Figure  4 

The  fluid  damping  ratio  Cp  plotted  against  V^;  from 
Griffin  & Koopmann  [19]. 


NO  SELF -EXCITED  OSCILLATIONS 


I 

I 


CASE  1:  Y»0  25.  Vr  =5.5  { ~ ^ | 

« » +204.7*  (W  *0.05) 
fi  * ARC  TAN  (Cdh/Cn,h^  *+201.8* 

Figure  5(a) 

Phase  distribution  of  the  flow-induced  forces  C . 

mh 

and  C , and  the  displacement  y;  measured  forces  from 
dn 

Sarpkaya  [21]. 
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+ COS  ar 


SELF- EXCITED  OSCILLATIONS 
POSSIBLE 


CASE2:Y.0.50,V,.5 
< = -26.5* 

ARC  TAN  (Cdh/Cmh)*-29.4* 

Figure  5(b) 

Phase  distribution  of  the  flow-induced  forces 

and  C.,  , and  the  displacement  y;  measured  forces  from 
dh 

Sarpkaya  [21].  * _ 
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Figure  6 

slncp  and  coae  plotted  agalnat  2Yj^  (scaled 
as  shovn) ; for  legend  see  Table  3. 
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Figure  7 

Cj^  sincpj^  and  -C^  sin€  plotted  against  2Yj^;  from 

Griffin  & Koopmann  [ 19]  - from  Sarpkaya  [21]  

from  Mercler  [20]  (all  circular  cylinders). 
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